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Our first-principles calculations demonstrate that Ceo-nB^ and Ceo-mNm molecules can be en- 
gineered as the acceptors and donors, respectively, which are needed for molecular electronics, by 
properly controlling the dopant number n and m in Ceo- As an example, we show that acceptor 
C48B12 and donor C48N12 are promising components for molecular rectifiers, carbon nanotube-based 
n-p-n (p-n-p) transistors and p-n junctions. 

PACS (numbers): 33.15.-e, 85.65. -l-h, 31.15.Ar 
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Modern microelectronics and computation are advanc- 
ing at an extremely fast rate because of remarkable cir- 
cuit miniaturization [1]. However, this trend will soon 
reach the scale of atoms or molecules. To continue toward 
faster and smaller computers, new schemes are required. 
Molecular electronics [2] is one such approach. 

One major problem in molecular electronics is con- 
necting the functional molecules and assigning the ob- 
served electrical properties in an unambiguous way to 
the molecules in question [3,4]. FuUerenes [5] are large 
enough to be identified by transmission electron mi- 
croscopy or scanning probe methods [4], are stable and 
easy to build into molecular circuits [5,6], and might be 
inserted into single-walled carbon nanotubes (SWNT) 
[7]. Hence, fuUerenes should be ideal components for 
molecular electronics. 

As in semiconductor electronics [1], acceptor/donor 
pairs are critical for use in molecular electronics, for ex- 
ample, molecular rectifiers [2], nanoscale p-n-p transis- 
tors and p-n junctions [8]. For traditional silicon doping, 
group V atoms (for example, phosphorous) act as donors 
and group III ( for example, boron) are acceptors. Anal- 
ogous acceptor /donor schemes are needed in molecular 
electronics [2]. FuUerenes are unique because they can 
be doped in several different ways (for example, endo- 
hedral [5], substitutional [5,9-11] and exohedral doping 
[5] ) . This should provide a wide range of possible accep- 
tor/donor schemes. 

To design active molecular devices, components with 
controllable electronic properties are needed [2,3] . For ex- 
ample, to obtain molecular rectification, the lowest unoc- 
cupied molecular orbital (LUMO) of the acceptor should 
lie at or slightly above the Fermi level of the electrode and 
above the highest occupied molecular orbital (HOMO) of 
the donor [2] . Hence it is important to search for desired 
acceptor/donor pairs which satisfy the requirement. In 
this letter, we suggest a controlled approach to obtain 
such pairs from Cgo molecules by using substitutional 
doping. Because the average carbon-carbon bond length 
in Ceo is slightly larger than that in graphite, which can 
only be substitutionally doped by boron, and the force 



constants [5] are somewhat weakened by the curvature of 
the Cgo surface, both boron and nitrogen can substitute 
for one or more carbons in Cgo [9-11]. Our first-principles 
calculations demonstrate that Cgo-nBn and Cgo-mNm 
molecules can be engineered as the acceptors and donors, 
respectively, which are desired for molecular electronics 
by properly controlling the dopant number n and m. As 
an example, we present the electronic properties of the 
acceptor C48B12 and donor C48N12, and discuss their po- 
tential applications in molecular electronics. 

First, we discuss the electronic properties of Cgo-nBn, 
Ceo-nNn and Cgo to show that Cgo-nBn and C6o-nN„ 
can act as controlled dopants. The number n ranges from 
1 to 12. Since there are many isomers [12] of C6o-nX„ at 
fixed n, we only take the dopant assignment of Hultman 
et al. [11] as an example [13]: each pentagon receives 
a maximum of one dopant X, and two X are separated 
from each other by two carbon atoms. The optimized 
geometry and total energy are calculated by using the 
Gaussian 98 program [14,15] with the B3LYP [16] hy- 
brid density functional theory (DFT) [17] and 6-31G(d) 
basis set. The calculated results are summarized in Ta- 
ble I. It is found that the binding energy of Cgo-nXn 
(X=B,N) decreases monotonically with increasing inte- 
ger n. For Ceo-nBn, the binding energy is 0.21 eV/atom 
to 0.38 eV/atom lower than that {Eb = 6.98 eV/atom) of 
Ceo- For Ceo-nN„, it is 0.10 eV/atom to 0.61 eV/atom 
lower than that of Ceo- The stabilities of Ceo-nB„ and 
Ceo-nN„ are comparable to but less than that of Ceo- 
In 1991, Smalley and coworkers [9] successfully synthe- 
sized Ceo-nB„ with 1 < n < 6. In 1995, Hummelen 
et al. [10] produced C59N. Very recently, Hultman et 
al. [11] have reported the existence of C48N12. These 
experiments show that certain stable Ceo-nX„ (X=B, 
N) can be made. Our results show that other B- or N- 
doped Ceo-nX„ structures for n < 12 have similar sta- 
bility. Our calculated ionization potential {Ej = 7.32 
eV) and electron affinity {Ea — 2.40 eV) for Ceo agree 
wcU with experiments (E/ = 7.54 ± 0.01 eV [18] and 
Ea = 2.689 ± 0.008 eV [19]), indicating the accuracy of 
our calculations for Ceo-nX„. Among Ceo and Ceo-nXn 
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(X=B,N) at fixed n, C6o-nB„ has the highest electron 
affinity, while Cgo-nNn has the lowest ionization poten- 
tial, indicating that Cgo-nBn and Cgo-nNn can serve as 
electron acceptor and donor, respectively. 

Calculations of LUMO and HOMO energies are nec- 
essary to explore combinations of donors and accep- 
tors suitable for molecular electronics. Table I shows 
the change of the LUMO/HOMO energies (and LUMO- 
HOMO gaps Eg = Ei — Eh) as we vary the dopant num- 
ber n in the Cgo-nXn molecule. Hence, acceptor/donor 
pairs which satisfy the required conditions in molecular 
electronics can be obtained by controlling the dopant 
number n and m in Ceo-nBn and Cgo-mNm. For ex- 
ample, based on the acceptor LUMO and donor HOMO 
shown in Table I, we may choose acceptor /donor pairs, 
such as C58B2/C58N2, C54B6/C54N6, C49Bn/C5iN9, or 
C48B12/C48N12 to build molecular rectifiers exhibiting 
similar rectification behaviors. In the following, we take 
C48B12/C48N12 as one example to show their applica- 
tions in molecular electronics. 

Table I: LUMO {Ei, in cV), HOMO {Eh, in eV), bind- 
ing energy {Eb, in eV/atom), ionization potential (Ei, 
in eV), and electron affinity {Ea, in eV) calculated for 
C6o-„X„ (X=B,N) using B3LYP/6-31G(d). 



n CgO-nBn 

El Eh Eh Ej Ea 

1 -4.26 -5.57 6.77 6.77 3.07 

2 -4.64 -5.25 6.75 6.35 3.53 

3 -4.21 -5.49 6.74 6.37 3.32 

4 -4.43 -5.60 6.73 6.79 3.32 

5 -4.59 -5.54 6.71 6.85 3.56 

6 -4.55 -5.49 6.70 7.01 2.98 

7 -4.49 -5.54 6.67 6.75 3.23 

8 -4.62 -5.34 6.66 6.44 3.50 

9 -4.65 -5.43 6.64 6.61 3.34 

10 -4.57 -5.44 6.63 6.57 3.45 

11 -4.45 -5.61 6.62 6.87 3.30 

12 -4.24 -5.58 6.60 6.73 3.08 



El Eh Eb El Ea 
-3.63 -4.72 6.88 5.93 2.15 
-3.71 -4.64 6.82 5.85 2.27 
-3.72 -4.69 6.78 5.92 2.29 
-3.91 -4.56 6.73 5.83 2.19 
-3.74 -4.56 6.69 5.95 2.20 
-3.74 -4.61 6.64 5.85 2.32 
-3.58 -4.59 6.60 5.91 2.39 
-3.60 -4.41 6.55 5.65 2.36 
-3.49 -4.45 6.51 5.77 1.90 
-3.31 -4.28 6.46 5.52 1.89 
-3.36 -4.35 6.41 5.71 2.29 
-2.61 -4.38 6.37 5.66 1.49 



The molecular geometry (size) and symmetry play im- 
portant roles in molecular electronics [2,3]. For example, 
an important factor in inducing rectification is some ge- 
ometric asymmetry in the molecular junction [3]. Cgo 
is a truncated icosahedron with a perfect Ih symme- 
try. Then, the changes in geometry (size) and symme- 
try of Cgo-nXn, due to the dopant-induced effects, are 
needed to know. Fig.l shows the optimized structure of 
C48B12. For comparison, similar calculations were done 
for Cgo and C48N12. The equilibrium C48B12, similar to 
C48N12 [11], has one dopant (boron) per pentagon and 
two dopants (boron) preferentially sit in a hexagon. In 
C48X12 (X=B, N), the doping- induced distortion from 
the perfect Cgo sphere is not localized to the neighbor- 
hood of each dopant but extends throughout the whole 



cage. The symmetry of C48B12 is the Ci point group, dif- 
ferent from the Sq symmetry for C48N12 [11]. This is due 
to their opposite electronic polarization, as discussed be- 
low. The molecular structure is examined by calculating 
the average distances (or radii R) from each atom to the 
center of the molecule. C48B12 has an ellipsoidal struc- 
ture, as in C48N12, with 10 unique sites labelled 1 to 10 in 
Fig.l and 10 different radii ranging from 0.347 nm (site 
2, 2') to 0.387 nm (site 8, 8'). In contrast, the 10 unique 
radii for C48N12 range from 0.340 nm (site 7) to 0.362 nm 
(site 5). For Cgo, each carbon atom has an equal radius 
R = 0.355 nm, the same as that found in experiment [20]. 
Furthermore, we find 15 unique bonds (specified in Fig.l) 
in C48B12: 6 boron-carbon bonds with lengths between 
0.154 nm and 0.159 nm, and 9 carbon-carbon bonds with 
0.139 nm to 0.150 nm lengths. In comparison, C48N12 has 
6 nitrogen-carbon bonds with lengths ranging from 0.141 
nm to 0.143 nm, and 9 carbon-carbon bonds with lengths 
from 0.139 nm to 0.145 nm. Cgo, however, has one kind 
of single C-C bond (0.145 nm) and one kind of double 
C=C bond (0.139 nm), which are in excellent agreement 
with experiments (0.14459 nm and 0.13997 nm [20]). 




FIG. 1.: C48B12 structure optimized with B3LYP/6- 
31G(d). Red (grey) is for B (C) sites. The 15 bonds 
between labeled vertices are all unique. 

To determine the electronic polarization, MuUiken 
charge analysis [21] was performed for the on-site charge 
Qm of atoms in C48B12, C48N12 and Ceo. C48B12 has 
two types of boron dopants: one boron with net Mul- 
liken charge Qm = 0.164 ]e] (site 5) and one with 
Qm = 0.187 \e\ (site 9), while carbon atoms in C48B12 
have negative Qm in the range from -0.004 \e\ (site 3) to 
-0.083 [e] (site 8). In C48N12, there are two types of ni- 
trogen dopants: one nitrogen with Qm = —0.595 \e\ (site 
5) and one with Qm = —0.600 \e\ (site 9), and two types 
of carbon atoms: one-fourth of the carbon atoms with 
negative Qm (-0.0125 [e] at site 2 and -0.0298 [e] at site 
3) and three-fourths of the carbon atoms with Qm in the 
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range of 0.192 |e| (site 6) to 0.227 |e| (site 4). In Ceo, 
Qni — f^or each carbon atom. Although the MuUiken 
analysis cannot estimate the atomic charges quantita- 
tively, their signs are correct [21]. Hence, Ceo is isotropic, 
while C48B12 and C48N12 have opposite electronic polar- 
izations, leading to their different symmetries. 
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FIG. 2.: Orbital energies and symmetries of the 
fcth eigenstate for C48X12 (X=B,C,N) calculated with 
B3LYP/6-31G(d). 

Applications of acceptors Ceo-nB^ and donors 
Ceo-nN„ require the proper line-up of energy levels. 
Upon boron or nitrogen doping, the degeneracy of 
the original Ceo levels, shown in Fig. 2, is removed by 
the structural distortion of the icosahedral symmetry. 
C48B12 (C48N12) is isoelectronic with C^jo^^ (Cgo^^)- Re- 
placing 12 carbon atoms by 12 boron atoms results in 
removing 10 electrons from the HOMO {jiu symmetry, 
Eh = —5.99 eV) of Ceo and 2 electrons from the 2nd 
HOMO (gg symmetry) of Ceo. In contrast, replacing 
12 carbon atoms by 12 nitrogen atoms leads to a com- 
plete filling of 6 electrons in the LUMO (ti„ symmetry, 
= -3.22 eV) of Ceo and 6 electrons in the 2nd LUMO 
{f,\g symmetry) of Ceo- Quantitatively, Fig. 2 shows the 
orbital energies of C48B12 and C48N12. They are differ- 
ent from each other because of the difference in valency 
of B and N atoms. For C48N12, the HOMO is a doubly- 
degenerate level of ag symmetry, while the LUMO is a 
non-degenerate level with au symmetry. For C48B12, the 
HOMO and LUMO are doubly-degenerate levels of au 
and ag symmetries, respectively. These calculations show 
that C48Bi2's LUMO is just above C48Ni2's HOMO and 
the approximate potential difference is about 1.63 eV. 
Thus, both C48B12 and C48N12 have precisely the prop- 
erties required for a proper rectifier behavior as will be 
demonstrated below. 

Various experiments [22] have shown that individual 
molecules can act as molecular rectifiers [2] . The inset in 
Fig. 3 (a) shows a rectifier molecule consisting of a donor 



C48N12 and a acceptor C48B12 connected by a tunnel- 
ing bridge (a cr-elcctron system, — (CH2)e— )• The bridge 
forms a potential barrier that isolates the donor from 
the acceptor on the time-scale of electron motion to or 
from the electrodes [2] (see review [6] about synthesis 
techniques which have successfully bridged Ceo dimers 
through electroactive spacers). Following the Aviram- 
Ratner scheme [2] , we expect that electron current would 
pass from the cathode to the anode via acceptor, bridge 
and donor. Since C48Bi2's LUMO is -4.2 eV, we may 
choose metallic carbon nanotubes (work function rj 5 
eV [23]), gold, silver or copper ( </) « 4.4,4.7,4.8 eV, 
respectively [24]) as the electrodes. As discussed be- 
fore, C48B12 and C48N12 have the properties required for 
proper rectifier behavior. An initial B3LYP/3-21G cal- 
culation for the full, connected fullerene/spacer/fuUerene 
system shows that the covalent bond between the spacer 
and doped fuUerenes shifts the electronic structure of 
the doped fullerenes. However, the HOMO and LUMO 
of the acccpter-spacer-donor system are still localized 
on the acceptor and donor sides, respectively, ensur- 
ing the desired rectifier behavior. The same properties 
hold for the acceptor/bridge/donor system [25]. Thus, 
a smaller threshold voltage for conduction is expected 
in one direction than in the other direction. As an ex- 
ample, we have investigated the rectifying characteris- 
tic of the C48Bi2/CeHi4/C48Ni2 system by using non- 
equilibrium Green function theory in conjunction with 
a density functional based tight-binding model [26,27]. 
The rectifier molecule is connected to gold contacts that 
we treat in the s-wave approximation with a constant 
density of states near to the Fermi level [28]. The calcu- 
lated current /applied potential characteristic is shown in 
Fig. 3 (a). A typical rectification characteristic is obtained 
with a turn-on bias close to the energy difference be- 
tween C48Bi2's LUMO and C48Ni2's HOMO. Recently, 
Joachim et al. [4] have successfully observed the electrical 
current flowing through an individual Ceo molecule with 
a scanning tunneling microscope (STM). In light of this 
progress and the calculated size of C48X12 (X=N,B), our 
rectifier molecule also should be identifiable by STM. 

Heterojunctions for molecular electronics can also be 
formed using Ceo-nX„ donors and acceptors. Hybrid 
nanostructures formed by filling a SWNT with Ceo [7] 
have been observed, and can be superconducting [7] or 
metallic [29]. Our DFT calculations, done within the 
local density approximation with double numerical ba- 
sis including d-polarization function [15,30], show that 
the incorporation of the C48B12 or C48N12 into a (10,10) 
[5] or (17,0) [5] SWNT is energetically favorable. About 
2.4 eV binding energy per molecule is gained after the 
C48B12 or C48N12 molecule is inserted periodically inside 
the SWNTs. In our calculations of doped SWNTs with 
a one-dimensional periodic boundary condition along the 
tube axis, three unit cells of a (10,10) tube or two unit 
cells of a (17,0) tube are included in the supercell with 
one C48B12 or C48N12. Charge analysis found that plac- 
ing an acceptor C48B12 into a (17,0) tube puts a -1-0.67 
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|e| charge on the SWNT, while incorporating a donor 
C48N12 into a (17,0) SWNT puts a -0.39 |e| charge on the 
SWNT. Similar results are obtained for a (10,10) SWNT. 
Hence, putting C48B12 into a semiconducting tube results 
in a p-type region on the SWNT, while filling of donors 
C48N12 into a semiconducting tube leads to a n-type re- 
gion on the SWNT. Thus, it is possible to use them to 
get n-p-n andp-n-p transistors [8]. As shown in Fig. 3(b), 
doping a (17,0) semiconducting SWNT with C48B12 and 
C48N12 should also make a SWNT-based p-n junction. 
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FIG. 3: (a) Calculated current through a rectifier 
molecule (inset) which consists of C48B12 (left, red for B 
atoms) and C48N12 (right, blue for N atoms) connected 
by C6H14 (middle, green for H atoms) and is coupled to 
two Au electrodes via S atoms on the two sides. The 
two C atoms in the ends of the C6H14 bridge are con- 
nected to the N (B) atoms in C48N12 (C48B12) respec- 
tively. A B3LYP/3-21G calculation shows that the con- 
nected rectifier molecular superstructure is stable relative 
to the separated, individual components by about 19.545 
kcal/mol energy, (b) a prototype for C48Xi2@(17,0) 
SWNT-based (X=B,N) p-n junction. 

In conclusion, we show that acceptor Cgo-nBn and 
donor Cgo-mNm pairs needed for molecular electronics 
can be obtained by properly controlling the number n 
and m of the substitutional dopants in Cgo. We demon- 
strated the rectifying characteristic of a molecular rec- 
tifier built from our engineered acceptor/donor pairs. 
Hetero junctions for molecular electronics can be made 
by inserting these dopants into semiconducting carbon 
nanotubes. Efficiently synthesizing those acceptor/donor 
pairs would be of great experimental interest within reach 
of today's technology. Very recently, small nitrogen- 
substitutionally-doped fuUerenes have been reported, [31] 
showing experimental progress in this direction. 
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